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The study and characterisation of transverse-momentum-dependent distribution functions 
(TMDs) is a major goal of the Electron-Ion Collider (EIC) physics programme. The study 
of gluon TMDs poses a greater challenge than for quark TMDs in DIS measurements, as 
gluons do not directly couple to photons. The study of D meson pairs has been proposed 
to provide access to gluon TMDs, but is demanding due to the rarity of D production. 
Here, we discuss the feasibility of such a measurement, and touch upon wider issues to be 
considered when measuring TMDs at the EIC. 



1 Introduction 

The EIC is proposed as a next-generation collider facility, with unprecedented luminosity and 
the ability to study both nucleons and nuclei at a variety of energies. There are two proposed 
realisations: eRHIC at Brookhaven National Lab (BNL), and ELIC at Jefferson Lab. Building 
on the legacy of HERA, RHIC, CERN and others, it will deliver a broad programme of nucleon 
and nucleon structure studies with extremely high precision. One of the key topics of interest 
at the EIC is the family of distribution functions known as TMDs. TMDs add a transverse mo- 
mentum component to conventional (one-dimensional) parton distribution functions (PDFs), 
providing a "3D" picture of partons in the the nucleon: TMDs essentially allow a tomographic 
imaging of the motion of partons at the femtoscale. Work at HERMES and COMPASS, among 
others, has provided tantalising hints about these distributions. However, data remains rela- 
tively scarce and a detailed understanding is far from achieved. The high luminosity and broad 
kinematic reach of the EIC will provide the perfect environment for not just exploring these 
functions, but characterising them in exquisite detail. 

The TMDs of quarks and gluons are of equal importance, but studying gluon distributions in 
electron-proton collisions is more challenging. Unlike electrically charged quarks, gluons do not 
directly couple to photons, which mediate the c-p interaction. A potential route to access gluon 
distributions is to study the production of pairs of D mesons. Here, we discuss the potential of 
such a method in measuring the gluon Sivers distribution. The Sivers distribution is a member 
of the TMD family that correlates the transverse momentum of partons in a nucleon with the 
transverse spin of the nucleon. It has enjoyed a great deal of theoretical and experimental study 
in recent years, and as such stands as an exemplar for the wider TMD family. 
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2 Studying the gluon Sivers function with D meson pairs 



The methodology for studying the gluon Sivers function with D mesons is detailed extensively 
in section 2.3 of pQ. For certain kinematics, a charm-anticharm pair produced from a gluon 
(via the process j*g —> qq) interacts like a single gluon. Thus, the QCD interactions that 
give rise to the Sivers distribution in gluons will give rise to the same Sivers distribution in 
the charm quark pairs. Quarks are not themselves detected, but their kinematics can be at 
least partly reconstructed from those of the D mesons they produce. As the Sivers distribution 
correlates spin and transverse momentum, a non-zero distribution results in a spin-dependent 
azimuthal asymmetry in the D pair transverse momentum distribution around the virtual 
photon mediating the DIS interaction. Measuring this asymmetry would allow the gluon Sivers 
distribution to be extracted. 

This signature is exciting as it provides an avenue to measuring the Sivers function where 
few (or no) others are currently known to exist. For example, dijet or dihadron production 
in proton-proton collisions, which would provide a natural environment for measuring gluon 
distributions due to the dominance of gluon interactions, are excluded due to the demonstration 
of factorisation breaking, precluding theoretical treatment. 

While exciting, the measurement is not without its difficulties. The rarity of D meson 
production, coupled with the necessity to detect and reconstruct not one but two in a single 
event, makes it statistically challenging. To investigate the feasibility of such a measurement, 
PYTHIA simulations were run to determine the event rates that could be expected at an EIC. 

3 PYTHIA simulation 

Charm events were generated with PYTHIA 6.416 for a high-energy EIC configuration: a 
20 GeV/c electron beam colliding with a 250 GeV/c proton beam. This configuration is most 
favourable for two reasons: an enhancement in the charm production cross section, and, at least 
in the BNL eRHIC design for the EIC, the maximum beam luminosity (see tabic [I]). 

The CT10 NLO PDF parameterisations 
[2] were used. Radiative corrections were not 
explicitly included. A cut on the DIS inelas- 
ticity of 0.01< y < 0.95 was imposed and 
a single Q 2 bin of 1 to 10 GcV 2 was used. 
Events were required to contain at least two 
D° mesons. The D° — > irK decay channel 
was used, due to the large branching ratio 
(3.87% [3]) and the relative ease of experi- 
mentally reconstructing it. Other D° decays 
are typically semileptonic or multi-hadon de- 
cays with three or more final-state products, 
significantly increasing the difficulty of, or ex- 
cluding, reconstruction. While other D mesons (D + , D~ , D *) could be used in an actual 
experiment, the increase in statistics from such decays was not sufficient to change the basic 
message of the Monte Carlo study: whether the analysis is feasible or not. Therefore, the 
D°— > ttK decay channel alone was used for simplicity. D°s generated directly from the DIS 
interaction and those generated via feed-down from heavier D mesons, for example D* + , were 
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Table 1: eRHIC design luminosities and charm 
cross sections (from PYTHIA) for different eR- 
HIC energy configurations (E eiectro „ x E proto „). 
Luminosity depends on the proton beam energy, 
but not the electron beam energy. 
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not distinguished. Studies of the kinematics of a D° compared to that of the heavier D generat- 
ing it found little difference. Thus, a feed-down D° is considered to provide as good a proxy for 
the charm quark kinematics as its parent. The correlation between charm quark and D meson 
kinematics was also examined. They were typically found to be closely correlated in azimuth, 
by A(j) <0.5 rad in the virtual photon frame. 

No detailed simulation of the effect of detector performance was done, but a simple ac- 
ceptance cut, requiring all particles to be more than one degree from the beam direction, was 
imposed. 



4 Results and discussion 



D° pairs were binned simultaneously in azimuthal angle around the virtual photon direction 
and transverse momentum, and the statistical uncertainty in each bin estimated as \/N (for 
N counts). Uncertainties were rescaled to represent 100 fb _1 of data. As the measurement is 
a single-spin asymmetry, error bars were also divided by proton beam polarisation, which was 
assumed to be 70% (the eRHIC design polarisation). The resulting statistical uncertainties are 
shown in Figure [l] for two bins in D meson pair transverse momentum and 10 bins in azimuthal 
angle. 

As seen in the figure, uncertainties of a 
few to several percent can be anticipated, 
depending on the pair transverse momen- 
tum range. The potential asymmetries de- 
rived from the model in [TJ depend on the 
the pair transverse momentum and can ex- 
ceed 20% for kp er p >1 GeV/c. This indi- 
cates that measuring D pair production is 
indeed feasible as an avenue to measuring 
the gluon Sivers function. At its maxi- 
mum luminosity, and conservatively assum- 
ing an operational efficiency of 50%, eRHIC 
would accumulate about 3 fb _1 per week. 
100 fb _1 would therefore require close to 
entire year of running at eRHIC. The eR- 
HIC physics programme would generally 
demand multiple energy configurations to 
be run in each year, so such an amount 
would likely represent a multi-year goal. 

The measurement of D mesons high- 
lights key performance requirements for an 
EIC detector that apply more generally to 
other semi-inclusive and exclusive measure- 
ments at an EIC. Perhaps most important 
is the requirement for precise tracking and 
particle identification (PID) over a wide 
range in rapidity. As is seen in Figure [2j 

not only are hadrons produced in the proton beam direction, but also at significant rates in the 
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Figure 1: Estimated statistical uncertainties on 
the spin asymmetry in D meson pair production 
with 100 fb _1 integrated luminosity for two bins 
in pair transverse momentum: 0.5 to 1 GeV/c 
(red, dashed) and 1 to 2 GeV/c (black, solid). Bin- 
to-bin variations are due to limited statistics in 
the Monte Carlo sample. 
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Figure 2: Distriubution of ir from D° decays in momentum and pseudorapidity (colour online). 



central region and the electron beam direction. Hadron momenta are typically modest in the 
central region (less than a few GeV/ c) , but extend to to several tens of GeV/ c in the proton beam 
direction. Therefore PID detectors capable of distinguishing hadrons up to high momentum 
will be required in this region. Furthermore, in common with inclusive events, the scattered 
electron should be detected and measured precisely to determine the DIS kinematics. An EIC 
detector should therefore provide high momentum and angular resolution and hadron-electron 
separation over a several units of rapidity. Whilst doing so, the material budget of the tracking 
detectors must be kept low enough to allow measurement of electrons down to low momenta 
(less than 1 GeV/c), where multiple scattering becomes an issue. 
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